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CHAPTER 1 

Problem Definition and Model Description 

The size of an ocean platform and the variation of 
environmental excitations on it make the determination of the 
integrity of the platform a vital concern to the ocean 
engineer. A method of integrity certification which is under 
development utilizes on-location measurement of the platform 
natural frequencies to determine if a structural failure has 
occurred. This technique requires an initial , accurate 
baseline determination of the platform natural frequencies. 
These baseline measurements will be compared to values 
obtained at later dates to establish if any changes have 
occurred. However, before a shift in natural frequency can 
be related to damage, all other sources of change, such as 
live platform load variations and marine growth, must be 
evaluated. In this paper, the author addresses the impact 
of a non-linear structural dynamic element on the measurement 
of the fundamental natural frequency. 

For example purposes, the platform to be analyzed is a 
deep ocean drilling platform with a centrally located, 
unattached, conductor which extends from the platform deck 
to the ocean floor. The purpose of the analysis is to 
determine the effect of the non-linear interaction between 
the platform and the conductor on the measured natural 
frequency of the entire system in the fundamental flexural 
mode. In order to limit the problem to a manageable size. 
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several constraints must be placed on the system. Principal; 
among these is the use of a periodic excitation to drive 
the system. This constraint simplifies the identification 
of the platform response frequencies. Secondly, the initial 
separation between the platform and the conductor is set to 
zero. In the author's opinion, this position presents a worst 
case condition, and it is also the most likely from engineering 
considerations of the conductor installation. The excitation 
loads were varied in order to demonstrate a breakpoint 
between the linear and non-linear behavior of the structure. 

The remaining system parameters were set equal to specific 
values based on engineering judgement and the best information 
available. The greatest uncertainty in parameter estimation 
occurred in the determination of the conductor parameters. 

A more definitive investigation of these quantities should 
be performed before further system analysis is conducted. 

The platform can be modeled as a mass, spring and damper 
system which surrounds another mass, spring and damper system 
modeling the conductor. The platform mass is free to move 
along its horizontal axis, and if sufficient excitation is 
provided, collide with the conductor' mass. This collision 
causes a reduction in the platform movement while increasing 
the motion of the conductor. The figure on the next page 
illustrates the model. 
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ml platform mass 

kl platform spring constant 

bl platform viscous damping coefficient 

m2 conductor mass 

k2 conductor spring constant 

b2 conductor viscous damping coeffidient 

xl platform displacement from equilibrium position 
x2 conductor displacement from equilibrium position 



FIGURE 1.1 SYSTEM SCHEMATIC DIAGRAM 
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A mathematical solution for the response of tht model is derived 
in the following manner; 

1. Equation of motion of the platform moving 
separatel y : 



miXi + bi^i + kjXi = F sin cot, 

or, XI = (l/mj) (F sin (ot - biXj - kjXj) (1) 

where x ^ is relative to the unforced equilibrixim position of the platform. 

2. Equation of motion of the conductor moving 
separately ; 

maX2 + b2X2 + k2X2 = 0 

or, X 2 = (l/m 2 ) (-b 2 X 2 - k2X2) ( 2 ) 

where x ^ is relative to the unforced equilibrium position of the conductor. 

3. Equation of motion for the combined platform 
and conductor masses : 



(mi + m 2 )xi + biXi + b 2 Xi + k^Xj + k 2 X, = F sin cot 
or, XI = (l/mi+ra 2 )(F sin cot - (bi+b 2 )xi - kixi - k 2 X 2 ) (5) 

In equat ions 1 , 2and 3, masses, damping constants and 
spring constants are denoted by m, b and k respectively. 
Subscript 1 designates the platform and subscript 2 denotes 
the conductor. Equations 1 and 2 completely describe the 
“»y g' p ^‘n-a-e* 'tfB f thc ‘ ~»y^-evt^hen- the- plat ferm stfhd-^ondy-ct-o^r— - 
in motion separately. Equation 3 describes the response 
after a collision has occurred. However, several questions 
about the system have yet to be answered. Specifically, 
what happens at the time of <ol1i.slon, how is (he t J me of 
collision determined, and how is the time of separation 
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determined? The answers to these questions provide the 
remaining equations which will define the motion of the 
platform/conductor system at all times. 

To deteimine when a collision occurs, the gaps between 
the platform and the conductor must be continuously determined. 
When the value of the gap distance on either side of the 
conductor goes to zero, a collision has occurred. The 
following equation defines this feature mathematically. 

4. Total gap distance; 
gapi + gaPa = do - di 

where do is the platform aperture diameter and di is the 
conductor outside diameter. 

5. Gap distance at any time, t: 

gapi(t) e gapj(t-l) + (XjCt) - XjCt-l)^ -(Xi(t) 

- xi(t-l)) (5a) 

gapaCt) = do - di - gapi(t) (5b) 

where (t-1) defines the time increment just prior to the 

time,t, being calculated. 

At the time of collision, energy will be lost. However, 
momentum is conserved. Therefore the velocity of the system 
can be determined as shown below. 
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6. Conservation of momentum: 
mix 1 = (m 1 + m 2 ) x 

or, X = (mi/mi + m 2 ) Xi (6) 

Since both the conductor and platform are moving at the same 
velocity^: the value of x is the value of Xi in equation 3. 

Finally, an equation is needed to determine if the plat- 
form and conductor have separated. The contact force can be 
used to perform this determination. The diagram below 
illustrates the contact force, f(t). 



f(+) 




m2 X 



FIGURE 1.2 CONTACT FORCE COMPONENTS 



7 . Contact force : 



f (t) = b2X + k2X2 



m2X 



(7) 
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In summary, the ocean platform has been modeled and the 
mathematical equations have been written to define the motion 
of the platform system at any point in time. The next 
chapter illustrates the process of non-dimensionalizing these 
equations which will lend them to analysis on a more general 
level. A numerical integration method of solution will then 
be presented after which the results of the analysis will 



be determined. 
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CHAPTER 2 

Non-dimensionalizat ion of the System Equations of Motion 

The removal of dimensions from the system equations of 
motion performs two important functions. First, the dimension- 
less parameters of the system have a greater application. 
Secondly, systematic variation of the dimensionless parameters 
will provide data which can be compared to known standards 
when the excitation is periodic. 

Eighteen parameters are contained in equations 1 through 
6 in chapter 1. The equations can be non-dimensional ized 
by proper selection of reference quantities for length, mass 
and time. The author defines F/ki, mi and 2 tt/woi as the 
reference quantities for length, mass and time respectively. 

The variables in equations 1 through 6 for displacement, 
velocity, acceleration, total gap distance and time are made 
dimensionless by division and differentiation using these 
reference values. The conductor parameters of mass, viscous 
damping coefficient and spring constant are non-dimensional- 
ized by dividing by the respective platform parameters. After 
substitution of the new variables into equations 1,2,3 and 
7, the damping coefficient of the platform is divided by 
product of two times mj times Wgj . This ratio is defined as 
the damping factor, , and is dimensionless. The contact 
force variable becomes dimensionless after all the parameters 
of equation 7 have been non - di men s i on n 1 5 z c 4 . 
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For clarity, the contact force variable is redefined F to 
identify it as dimensionless. The remaining parameter, 
excitation frequency, is nondimensionalized by dividing by 
the platform natural frequency. The dimensionless parameters 
are summarized in Table 2.1. 



Table 2.1, Summary of Non-Dimensional Parameters 



Non-dimensional 


Dimensional 


Variable 


Ratio 


r 


Xi/(F/ki) 


c 


X2/(F/ki) 


T 


t/ (2 tt/ojo 1 ) * 


P 


kz/ki 


y 


m 2 / j 




b 2 /b 1 




bi/2miWoi 


Yi 


W/Uo 1 


d| 


l/(F/ki) dxi 


dr 


l/(F/ki) d^xi 


d5 


l/(F/ki) dX2 


d c 


l/(F/ki) d^X2 


dx 


1/(2it/(Ooi) dt 


dt ^ 


1/ (4ir^ /w? 1 ) dt^ 


f 


F/mi g 


D 


(do - di)/(F/ki) 


*Wo 1 = (ki/mi) ^ 





The equations of motion from chapter 1 and their dimensionless 
equivalents are presented on the following pages. 
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1. Platform motion moving separately: 

Xi = (l/mj) (F sin cot - bjXj - kjXj) (1) 

NON- DIMENSIONAL I ZED 

f = 4ir^(sin Zttyit- r) - AuCif (1*1) 

2. Conductor motion moving separately: 

xj = (l/mz) (-b2S:2 - ^ 2 X 2 ) (2) 

NON-DIMENSIONALIZED 

c = (1/y) (4Tn{;5ii - 4 tt^pc) (2.1) 

3. Combined platform/ conductor motion: 

Xi = (l/mi+m 2 )(F sin cot -(bi-*b 2 )xi - kiXi-k 2 X 2 ) (3) 
NON-DIMENSIONALIZED 



r = (1/1+y) (4 tt^ (sin 2 ttyiT - r - pc) -47T(Cif+ii»Cif)) (3.1) 

4. Gap distances at any time,t: 

gapi(t) = gapi(t-l) + (X 2 (t) - X 2 (t-]) - 



(xi(t) - xi(t-])) (5a) 

gap2(t) = do - di - gapi(t) (5b) 

NON-DIMENSIONALIZED 

gapi ( t) = gapi(T-l) + (c(t) - c(t- 1 ))- 

(r(x) - r(T-l)) ■ (5a.l) 

gap 2 (t) = D - gap 1 (t) _ 

5. Conservation of momentum: 

X = (mi/mi + mi) Xj (6) 

NON-DIMENSIONALIZED 

r = (1/(1 ) r (6.1) 
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Equation 6,1 appears inconsistent. However, it should be 
interpretted as the velocity the moment after a collision is 
equal to a constant times the value of the velocity just 
before the collision. 



Contact Force 




f(t) = b2X + kjXj - nijX 


(7) 


Non-dimensional ized 




f = 'KiT + pc - (l/4ir^)vir 


(7.1) 



The equations of motion describing the platform/conductor 
system have been non-dimensional ized and are in a form which 
lends then to analysis through the use of an integrative 
method of solution. 
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CHAPTER 3 

Method of Solution and Error Analysis 

Since the interaction between the platform and the 
conductor are quite non-linear, analytic methods of solution 
must be bypassed for numerical means. Of these, the Runga- 
Kutta (R-K) is most fitting for this analysis because the 
method is se 1 f - st art ing and easily adaptable to computer 
programming . 

Classically, the R-K method expresses the difference 

between the values of y at x(n+l) and x(n), where x(n+l) 

andx(n) are (n+1) and (n) increments of x away from a reference 

point Xq, if the increments are constant as they will be for 

this analysis. Mathematically this statement becomes: 

m 

yCn+1) - y(n) = Ea^k^ (8) 

where the a.'s are constant, the k 's are determined by Taylor 
1 1 

series expansion and 'm' is the order of the R-K method being 
used. Cursory examination of equation 8 indicates that higher 
orders of the R-K method yield a greater number of terms for 
each value of y(n + l) . It is to be shown that h'.gher orders also 
yield greater accuracy. 

The principle error in the Runga-Kutta method is the 
truncation error which is defined: 

T = c(h)"'’^^ + o(h)'"'"^ 

ID 

where: ra is the order of the R-K method 
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h is the increment (x^ - Xg ) 

c is a constant of f/dx™*^ . 

The order of the term is generally smaller than the h™'*'^ 

term and will not be considered further. The truncation error 

for any term may be defined: 

T^ oc • (number of steps) , 

where (number of steps) = constant/h. 

Therefore , 

T. « h*" 

1 

An estimation of the magnitude of the truncation error for a 

term can be determined by comparison of the integration values 

using spacing increments of h and 2h. 

T ^ y(x.) - y^^^^= c(2h)"’ = c(h)"*2"* 
m 11 

T = y(x.) - y^^^ = c(h)’", 
m 1 i 

where (2h) and (h) distinguish between integration increment 
values. Subtraction of the second equation from the first 
yields : 

y (h) . y.C2h) ^ c(h)’"(2"’ - 1) or, 
i 1 

T = (y Ch)_ y (2h) - 1) . 

mi i 

t> y. 

After actual calculation, T = 0.0.... b, where b is the m 

m 

decimal place. 

Keying on minimizing truncation error, the author selected 
a fourth order Runga-Kutta method for use in this paper. 

The numerical approximation equations for this method are as 
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fol lows : 

y(n + l) - y(n) = 5 (ki + 2k2 + 2ks + ki,) 

ki * h f (x(n) , y(n) ) 

kj = h f(x(n) + hh, y(n) + hk^) 

kj = h f(x(n) + hh, y(n) + ^kj) 

ki, = h f(x(n) + h, y(n) + kj) 

To perform the R-K integration for the determination of 
the platform and conductor displacements, it is necessary to 
transform the second order differential equations of chapter 2 
into first order differential equations. This elementary 
procedure yields the equations shown below. 



Platform Displacement (moving separately) 




Tj * 4 tt^( sin 2tt t - tj) - 4irCir2 



Conductor Displacement (moving separately) 




Cj = (1/p) (4Tnf/C,C2 - 47T^pCj) 




( 1 . 2 ) 

(1.3) 



( 2 . 2 ) 

(2.3) 



H)> 



Combined Platform/Conductor Displacement 



^2 = (1/1+y) (4TT^(sin 2 ttYiT - Tj - pCj) 

-4irCCir2 + 

= -a 

Platform/Conductor Gap 

gapj (t) = gapj(T-l) + (Cj(T) - Cj(T-l)) 
- (1*1 (t) - Tj (t-1) ) 

Co ntact Force 

= + pcj - Cl/4TT^)yr2 



(3.2) 

(3.3) 



(5a. 2) 
(7.2) 



The computer program utilized in this paper was designed 
and written by the author. A description, user's guide and a 
printout of the program with output are presented in Appendix 
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CHAPTER 4 
Analysis Results 

As it was mentioned in chapter 1, the problem has 
been constrained to maintain manageability. As the first 
iteration toward understanding the response characteristics 
of the platform system, the author set the dimensionless 
ratios of the conductor to platform mass, viscous damping 
coefficient and spring constant to constant values. The 
damping factor, Cj, and the total gap distance were also 
made equal to constants. The excitation force parameter, FOl, 
was given three discrete values for the purpose described in 
chapter 1. The ratio of excitation frequency to platform 
natural frequency, Yi» was given a range of seven values from 
0.5 to 5,0. This range yields a summary plot of the platform 
displacement which is similar in appearance to a standard 
transfer function plot. The values of the platform mass 
and natural frequency were estimated from previous platform 
analysis experience. The damping factor was also determined 
in this manner. With the values for natural frequency and 
mass, the platform spring constant was calculated. The ratios 
between conductor and platform parameters were then estimated. 
Using these percentages, the values of conductor mass, 
viscous damping coefficient and spring constant were calculated. 
The values of the system parameters are listed in Table 4.1. 



TABLE 4.1, System Parameter Values 



Parameter 


Value 




0.009998 


C 


0.03 


y 


0.1 


Y: 


0.5, 0.8, 1.0, 1.5, 2. 0,3. 0,5.0 


P 


0.01 


FOl 


0.01,0.001,0.0001 




2.07 X lO"* (lb-sec^) /inch 


m 


2.07 X 10® (lb-sec^) /inch 


bi 


3.12 X 10“ (lb-sec)/ft 


b, 


3.124 X 10® (lb-sec)/ft 




9.08 X 10“ Ib/in 


k2 


9.08 X 10^ lb/ in 




2.09 radians/second 




0.628 radians/second 


do 


19.5 inches 


di 


17.2 inches 



The system was analyzed for the platform response over the 
range of Yi for each excitation force. 

The output of the computer program was two plots. 

Each plot was coded on the abscissa to denote the type of 
response that was plotted and the specific values of Yi and 
FOl used for that computer run. A complete description of the 
coding and the individual displacement vs. t plots for the 
twenty-one parameter combinations are presented with their 
analyses in Appendix B. Each analysis discusses the effect 
of the excitation frequency on the response frequency, the 
dynamic equilibrium position either positive or negative of 
the static equilibrium zero position, the steady-state 
magnitude of the response and the primary response frequen- 
cies. 



As an example, if Yj and FOl were set to values of 
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0.5 and 0.001 respectively, the resulting platform displace- 
ment and fourier transform are shown in Figures 4.0 and 4.01. 
The platform response has a period of 2 t which is the period 
of the excitation frequency, 1 radian/ second . The dynamic 
equilibrium position is slightly negative. The Fourier 
transform yields two frequencies, the excitation frequency 
and the platform natural frequency. 




Figure 4.1, Platform Displacement V5. Yi 
In Figure 4.1 the dimensionless platform displacement 
is plotted against the ratio of excitation frequency to 
platform natural frequency, v,. The non-linearity of 
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system is demonstrated by the reduction in magnitude as the 
excitation force is decreased from 0.01 to 0,001, The 
decrease has also produced a frequency shift in the platform 
response to the right. The coincidental displacement plots for 
the excitation force equal to 0.001 and 0.0001 imply that 
the system is linear for these load levels. Thus, there is a 
breakpoint in the excitation force level below which the 
non-linear influence of the conductor on the platform natural 
frequency is negligible. However, the non-linearity has 
reduced the magnitude of the response. 

On the following pages. Figures 4.2 through 4.8 present 
plots of the absolute value of the Fourier transform of the 
platform displacement against the ratio of response frequency 
to the platform natural frequency which has been designated 
Y*. Each figure displays the displacement magnitudes for a 
specific value of the ratio of excitation frequency to platform 
natural frequency, y . The congruency of the plots for FOl 
equal to 0.001 and 0.0001 over all values of y^ reinforces 
the statements made previously concerning the non-linear 
effect of the conductor. 
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Figure 4.2, I R(w) I vs. Y*; Yi=0*5 
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Figure 4.3, |R(aj)I vs. Y*J Yi‘‘ 0*8 
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Figure 4 . 4 , |r((o)I vs. y*; Yi= 1*0 
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Figure 4.5, 1r(w)| vs. y *; Yx= 1*5 
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Figure 4.6, |R(oj)| vs. y*; Yi= 2.0 
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Figure 4.7, |R(w)| vs. y*; Yi= 3.0 
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Figure 4.8, 1r(cj)1 vs. Y*; Yi= 5.0 
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Principally, the analysis has shown that for large 
values of excitation force the system is non-linear. This 
is clearly displayed by the shift in the platform natural 
frequency. if the excitation load is below an as yet undeter- 
mined level, the natural frequency of the platform is unaffect- 
ed by the non-linear influence of the conductor. However, 
the non-linearity reduces the steady-state platform magnitude. 
The analysis consistently attributes the response frequency 
of the platform to a function of excitation frequency and 
the platform natural frequency. Normally, the platform 
would respond at the excitation frequency in the steady- 
state mode. However, the analysis was conducted by numerical 
integration until steady-state was attained. Therefore, when 
the Fourier transform of the platform displacement was taken, 
the transient effects of the problem "start-up” were included 
and presented in the transform results. This transform is 
more desirable than a pure steady-state transform since 
actual measurements on the platform are taken during random 
excitat ionand the transient effects will not have an oppor- 
tunity to die out. 

With the dimensionless excitation force equal to G.Ol, 
the dynamic equilibrium of the system shifted markedly 
negative if the excitation frequency was near the platform 
natural frequency. At frequencies above the platform natural 
frequency, the dynamic equilibrium always shifted positively, 
but the shift was nearly indistinguishable. The dynamic 
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equilibrium position^shift was inconsistent when' the exci- 
tation frequency equalled the platform natural frequency. 

The period of the maximum platform displacement in steady- 
state was coincidental with the excitation period throughout 
the analysis. No further conclusions can be drawn from the 
data . 

The Runga-Kutta numerical integration method has been 
an effective means of determining the dynamic characteristics 
of the system. It is evident that this analysis has raised 
more questions than it has answered because of the imposed 
parameter constraints. Even with these limitations, if several 
more of the parameters, i.e., the platform damping factor 
and the conductor/platform ratios were varied, the system 
response characteristics could be more accurately determined. 
The system should then be analyzed using an excitation approx- 
imating the environmental conditions to which the platform will 
be subjected to determine the fundamental flexural baseline 
response . 

The computer program has some minor problems which 
could be eliminated for greater efficiency. The program 
was designed to analyze the system with one set of parameter 
values per computer run. However, it would be more economical 
to have the capability to increment those parameters chosen 
for variation over a set of values. This improvement would 
keep the program on the computer until all of the variations 
had been completed. The Runga-Kutta section of the program. 
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although correct and complete, does not allow for varying 
the incrementing variable. If the difference between two 
adjacent, integrated values of a variable reaches a magnitude 
that is too large, the program may overflow the computer 
memory capacity and generate output which contains only error 
messages. The author endured this perennial set-back 
because of the input requirements of the Fourier transform 
and plotting subroutines. 
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The computer program was written to utilize a 4^^ order 
Runga-Kutta (R-K) numerical integration approximation for 
determining the dimensionless displacement of an ocean 
platform with a centrally located conductor. The values of the 
displacement are then printed out, plotted as a function of T 
and subjected to a Fourier transform subroutine whose output 
absolute values are plotted as a function of the to determined 
in the transform. Because the equations of the program are 
non-dimensionalized, the values of w used for the Fourier 
transform plot are a function of x, not real time, t. To 
convert w(fourier transform) to an co(real time) , one must 
divide the product of w(fourier transform) times toCplatform 
natural frequency) by 2 tt . The program has been designed to 
compute the platform response for one set of input parameter 
values per program execution. Thus, if Yi were assigned a 
range of values, the user must run the program on the computer 
once for each value of Yi* 

The program is initiated by setting the arrays for the 
platform displacement and velocity, the conductor displacement 
and velocity, the gap values and time equal to zero. This 
ste) is folloved by individually setting the values for the 
dimensionless parameters. The values of the dimensional 
parameters for excitation force, platform spring constant, 
platform aperture diameter and the conductor outside diameter 
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are then entered into the computer to calculate the value 
of the dimensionless total gap distance. The final initial 
conditions are input when the individual gap distances are 
set and the program begins execution. 

The analysis begins by computing the first integration 
values of platform displacement and velocity, and the 
conductor displacement and velocity using the 'masses 
separated’ R-K integration equations. Gapj is then calculated 
by subtracting the net change in platform displacement from the 
net change in the conductor displacement and adding this result 
to the old value of gapj. This new gapj is subtracted from 
the total gap distance to determine the new gap2. Gap^ and 
gapj are checked for a zero or negative value indicating a 
collision between the platform and the conductor. If no 
collision has occurred, time is incremented and the integration 
takes place again. If there has been a collision, the value 
of the conductor velocity is changed to equal a new value of 
the platform velocity which is determined by conservation of 
momentum. The contact force between the platform and the 
conductor is checked for a zero or negative value. If the 
contact force is zero or negative, the program returns to the 
'masses separated' R-K integration. If the contact force is 
greater than zero, the program executes the R-K approximation 
using the equations of motion for the combined masses. The 
program again checks for a zero contnct force and increments 
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the R-K approximations for separated or combined masses as 
the situation dictates until the limit of the increment count 
is reached. 

At this time the platform displacement is transferred 
into a plotting subroutine for a pictorial display of the 
result. The platform displacement array is then transformed 
into a complex number array and transferred to the Fourier 
transform subrout ine ( FFT) . The output of the FFT is transformed 
into an array of absolute magnitudes and returned to the 
plotting subroutine for display. Upon completion of the 
plot, the program ends. 



/ 
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USER'S GUIDE 

NOTE: This computer program has been designed for execution 

on the INTERDATA M70/M80 computers of the M.I.T. Joint 
Mechanical/Civil Engineering Computer Facility exclusively. 
Alterations will be required to fulfill the conditions of 
plotting and Fourier transform subroutines at other facilities 



The computer program is virtually self-contained requir- 
ing only two cards for execution per computer run once 
the parameter value cards have been installed. These cards 
provide the nomenclature for the 'x' and 'y' axes of the 
plotting subroutine. These cards are placed into the card 
deck as shown below in the following order: platform displace 

ment vs. t; absolute value of platform displacement vs. co . 

The format of these two data cards is 2A40. 

The platform/conductor system parameters are included in 
the program main body. Each parameter may be varied individ- 
ually by changing the appropriate card value. 

The program execution is controlled by submitting the 

card deck as follows: // JOB 

// FOR 

PROGRAM DECK (including subroutines 
RETURN 
END 

// XEQ2 
*EL FFTS 
*BC 2866 

DATA CARDS (TWO) 

// END 

The following pages list the program and an example 



output . 
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COMPLEX PWRC512) ,WTAB(256) 

INTEGER* 2 XLABC40) 

EXTERNAL AREA 

DIMENSION GAP1BR(512) ,GAP2BR(512) ,A(2,512), XSCL(4) 

COMMON Rl(512) ,R2(512) ,C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GAMMA1,PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2K0,C2K1,C2K2,C2K3,L.DAMC01,DELTAU,F01 
C SET INITIAL VALUES FOR ARRAYS AND PARAMETERS 

DO 1 1= 1,512 
R1(I) = 0.0 
R2(I) = 0.0 
C1(I) = 0.0 
C2(I) = 0.0 
GAPIBR(I) 5 0.0 
GAP2BR(I) = 0.0 
A(1,I) = 0.0 
A(2,I) - 0.0 
1 TAU(I) = 0.0 
DAMCOl =0.03 
L = 0 

PI = 3.14159625 
GAMMAl = 1.0 
DO = 19.5 
DELTAU = .1 
FOl = .01 



FO = 79984. 
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AKl = 90784. 

DI = 17.2 
RHO = .01 
PSI = .0099998 
AMU = .1 
J = 0 

DEAR = (DO - DI)/(F0/AK1) 

GAPlBR(l) = 0.0 
GAP2BR(1) = DEAR 
WRITE (5, 30) 

30 ^ORMAT(T10, 'TIME' ,T24,' PLATFORM DISP' ,T43, 'PLATFORM VEL' ,T60, 'COND 
2UCTOR DISP' ,T80,' CONDUCTOR VEL' ,T99, 'GAPl BAR',//) 

C START THE PROBLEM BY INTEGRATING WITH A 4TH ORDER RUNGA-KUTTA APPR 

N = 1 

4 DO-3' I U N,511 
J = I + 1 

TAU(J) = TAU(I) + DELTAU 
A(2,J) = TAU(J) 

CALL RKO 
CALL RKl 
CALL RK2 
CALL RK3 

R1(J) = AREA (R1(I),DELTAU,R1K0,R1K1,R1K2,R1K3) 

A(1,J) = R1(J) 

R2(J) = AREA ( R2 ( I) , DELTAU, R2K0, R2K1 ,R2K2 ,R2K3) 
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C1(J) = AREA (C1(I),DELTAU,C1K0,C1K1,C1K2,C1K3) 

C2(J) = AREA ( C2(I),DELTAU,C2K0,C2K1,C2K2,C2K3) 

C CHECK THE GAP DISTANCES BETWEEN THE PLATFORM AND THE CONDUCTOR 

GAPIBR(J) = GAPIBR(I) +(C1(J) - C1(I)) - (RU'J) - R1(I)) 

GAP2BR(J) = DBAR - GAPIBR(J) 

IF(GAP1BR(J) .LT. 0.0) GO TO 600 
700 IF(GAP1BR(J) .EQ. 0.0) GO TO 100 
IF(GAP2BR(J) .LT. 0.0) GO TO 800 
900 IF(GAP2BR(J) .EQ. 0.0) GO TO 100 
3 CONTINUE 

100 IF( J .EQ. 512) GO TO 400 

C WHEN A GAP GOES TO 0.0, CONSERVE MOMENTUM 

N = L + 1 
WRITE(5,40) TAU(N) 

40 FORMAT (/,T6, 'CONTACT OCCURRED AT TIME, ' ,2X,1PE12.4,/) 

IF(GAP1BR(N) .EQ. 0.0) GAP2BR(N) = DBAR 
IF(GAP2BR(N) .EQ. 0.0) GAPIBR(N) = DBAR 
R2(N) =(!./(!. +AMU)) * R2(N) 

C2(N) = R2(N) 

C INTEGRATE THE SYSTEM WITH COMBINED MASSES 

DO 5 1= N,511 

C CHECK THE CONTACT FORCE BEFORE INTEGRATION TO ENSURE NO SEPARATION 

FI = 2.* PI * GAMMAl * TAU(I) 

C ACCEL IS A DUMMY VARIABLE FOR THE ACCELERATION AT ANY TIME TAU 

ACCEL = (4. * PI/(1. + AMU)) * (PI * SIN(FI) -DAMCOl * R2(I) - PS 
21 * DAMCOl * R2(I) - PI * RHO * C1(I)) 
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C FHAT IS THE VARIABLE FOR THE CONTACT FORCE 

FHAT = FOl * PSI * DAMCOl * R2CI)/PI + RHO * FOl * Cl (I) - AMU * 
2F01 * ACCEL/ (4. * PI**2) 

IF( FHAT .EQ. 0.0) GO TO 15 
IF( FHAT .LT. 0.0) GO TO 15 
L = I 
J = L + 1 

TAU(J) = TAU(I) + DELTAU 
A(2,J) = TAU(J) 

CALL RKOl 
CALL RLll 
CALL RK21 
CALL RK31 

R1(J) = AREA (R1(I) , DELTAU, R1K0,R1K1,R1K2,R1K3) 

A(1,J) = R1(J) 

R2(J) = AREA (R2(I) , DELTAU, R2K0,R2K1.R2K2,R2K3) 

CICJ) = Cl(I) +(R1CJ) - RlCI)) 

C2CJ) = R2(J) 

GAPIBR(J) = GAPIBR(I) 

GAP2BRCJ) = GAP2BRCI) 

A 

GO^O 20 

15 WRITEC5,25) TAU(I) ,R1(I) ,R2CI) ,C1(I) ,C2(I) ,GAP1BR(I) 

25 F0RMAT(T6,1PE12.4,T25,1PE12.4,T43,1PE12.4,T61,1F'E12.4,T81,1FE12.4, 
2T97,1PE12.4) 

IF(J .EQ. 512) GO TO 400 
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WRITE (5, 11) TAU(J) 

11 FORMAT (/,T6, 'THE CONTACT FORCE IS ZERO;MASS SEPARATION OCCURRED 
2AT TIME, '2X.1PE12.4,/) 

N = L + 1 
GO TO 4 

20 WRITE(5,10) TAU(J),R1(J),R2(J),C1(J),C2(J),GAP1BR(J) 

10 F0RMAT(T6,1PE12.4,T25,1PE12.4,T43,1PE12.4,T61,1PE12.4,T81,1PE12.4, 
2T97.1PE12.4) 

5 CONTINUE 
400 M = 0 
GO TO 53 

54 DO 50 I = 1,512 

50 PWR(I) = CMPLX(R1(I) ,0.0) 

NBITS = 9 

NEWTB = 1 
INV = 1 

CALL FFT CPWR, NBITS, INV, WTAB, NEWTB) 

C THE VARIABLE 'S' IS THE INCREMENTAL FREQUENCY OVER WHICH THE 

C FOURIER TRANSFORM WILL BE PLOTTED, AND THE CONSTANT IN THE 'S' 

C EQUATION IS EQUAL TO ' 2*PI/ (512 . *DELTAU) ' 

S = 0.0 

DO 51 I = 1,512 ^ 

A(1,I) = CABS(PWR(I)) 

A(2,I) = S 

51 S = S + .6136 



53 NROW = 2 
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XSCL(l) =0.0 
XSCL(2) = 0.0 
XSCL(3) =0.0 
XSCL(4) =0.0 
NVARS = 1 
NPTS = 512 
NX = 2 
MOVE = 0 
LABEL = 4 
ISCL = 1 
FTIME = 0.0 
LOOK = 1 

READ(8.45) (XLAB(I), I = 1,40) 

45 FORMAT (40A2) 

CALL P ICTR (A , NROW , XLAB , XSC L , NVARS , NPTS , NX , MOVE .LABEL, ISCL, FTIME, 
2 LOOK) 

M = M + 1 

IF(M .EQ. 1) GO TO 54 
GO TO 35 

600 GAPIBR(J) =0.0 
GAP2BR(J) = DEAR 
C1(J) = R1(J) 

GO TO 700 

800 GAP2BR(J) =0.0 
GAPIBR(J) = DEAR 
C1(J) = R1(J) + DEAR 



GO TO 900 
35 STOP 
END 
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C THE RKO SUBROUTINE DETERMINES THE K1 COEFFICIENTS DURING 

C COMPONENT SEPARATION 

SUBROUTINE RKO 

COMMON R1(512),R2C512) ,C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GAMM.U,PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2KO , C2K1 ,C2K2 , C2K3 , L, DAMCOl , FOl , DELTAU 
PI2 * PI **2 

FI = 2.*PI * GAMMAl * TAU(L) 

RIKO = R2(L) 

R2KO = (4. *PI2) * CSIN(FI) - R1(L)) - 4. * PI * DAMCOl * R2(L) 
CIKO = C2(L) 

C2KO = -(4. * PI * PSI * DAMCOl * C2(L)/AMU) - (PI2 * 4. * RHO * 
2C1(L)/AMU) 

RETURN 



END 
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C THE RKl SUBROUTINE DETERMINES THE K2 COEFFICIENTS DURING 

C COMPONENT SEPARATION 

SUBROUTINE RKl 

COMMON R1(512),R2(512),C1(512),C2(512),TAU(512),PI AMU,GAMMA1 ,PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2 KO , C2 K1 , C2K2 , C2 K3 , L , DAMCOl , FOl , DELTAU 
PI2 = PI **2 
B = DELTAU / 2. 

T = TAU(L) + B 

El = R1(L) + B * RIKO 

E2 = R2(L) + B * R2K0 

E3 = C1(L) + B * CIKO 

E4 = C2(L) + B * C2K0 

FI = 2. * PI * GAMMAl * T 
RlKl = E2 

R2K1 = (4. * PI2) * (SIN(FI) - El) - 4. * PI * DAMCOl * E2 
ClKl = E4 

C2K1 = -(4. * PI*PSI*DAMC01*E4/AMU) - (PI2 * 4 . *RHO*E3/AMU) 

RETURN 



END 
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C THE RK2 SUBROUTINE DETERMINES THE K3 COEFFICIENTS DURING 

C COMPONENT SEPARATION 

SUBROUTINE RK2 

COMMON R1C512) ,R2(512) ,C1(512) ,C2(512) ,TAU(512) ,PI ,AMU,GAMMA1 ,PSI , 
2RH0,R1K0,R1K1,R1K2,R1K3,C1K0,C1K1,C1K2,C1K3,R2K0,R2K1,R2K2,R2K3, 
3C2K0,C2K1,C2K2,C2K3,L,DAMC01,DELTAU,F01 
B e DELTAU/2. 

PI2 = PI **2 
T = TAU(L) + B 
El = R1(L) + B * RlKl 

E2 = R2(L) + B * R2K1 

E3 = C1(L) + B * cm 

E4 = C2(L) + B * C2K1 

FI = 2. * PI * GAMMAl * T 
R1K2 = E2 

R2K2 = C4. * PI2) * CSIN(FI) - El) - 4. * PI * DAMC01*E2 
C1K2 = E4 

C2K2 = -(4. * PI * PSI * DAMCOl *E 4/AMU) - (PI2*4 . *RHO*E3/AMU) 
RETURN 



END 
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C niE RK3 SUBROUTINE UETEKMINES THE K4 COEFFICIENTS DURING 

C COMPONENT SEPARATION 

SUBROUTINE RK3 

COMMON R1(512),R2(512),C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GAMMA1,PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2K0,C2Kl,C2K2,C2K3,L,DAMCOl,DELTAU,F01 
B = DELTAU 
PI2 = PI **2 
T = TAU(L) + B 
El = R1(L) + B * R1K2 

E2 = R2(L) + B * R2K2 

E3 = C1(L) + B * C1K2 

E4 = C2(L) + B * C2K2 

FI = 2. * PI * GAMMAl * T 
R1K3 = E2 

R2K3 = (4. * PI2) * (SIN(FI) - El) - 4 . *PI*DAMC01*E2 
C1K3 = E4 

C2K3 = -(4.*PI*PSI*DAMC01*E4/AMU) - (PI2* 4 . *RHO*E3/AMU) 

RETURN 



END 



A-16 



C THE RKOl SUBROUTINE DETERMINES THE K1 COEFFICIENTS FOR THE 

C COMBINED MASSES 

SUBROUTINE RKOl 

COMMON Rl(512) ,R2(512) ,C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GAMMA1,PSI, 
2RH0,R1K0,R1K1,R1K2,R1K3,C1K0,C1K1,C1K2,C1K3,R2K0,R2K1,R2K2,R2K3, 
3C2 KO , C2 K1 , C2 K2 , C2 K3 , L , DAMCOl , DE LTAU , FOl 
FI = 2. * PI * GAMMAl * TAU(L) 

PI2 = PI **2 
RIKO = R2(L) 

R2KO = C4.*PI/(1.+AMU)) * (PI*SIN(FI) -DAMC01*R2 (L) -PSI*DAMC01*R2 
2(L) - PI * RHO * C1(L)) 

CIKO = C1(L) 

RETURN 



END 



A-17 

THE RKll SUBROUTINE DETERMINES THE K2 COEFFICIENTS FOR THE 
COMBINED MASSES 
SUBROUTINE RKll 

COMMON R1C512),R2(512),C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GAMMA1,PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2K0,C2Kl,C2K2,C2K3,L,DAMCOl,DELTAU,F01 
B = DELTAU + B 

T = TAU(L) + B 

FI = 2. * PI * GAMMAl * T 
El = R1(L) + B * RIKO 

E2 = R2(L) + B * R2KO 

E3 = C1(L) + B * CIKO 

RlKl = E2 

R2K1 = (4. * PI/C1.+ AMU))*(PI*SIN(FI) -DAMC01*E2/PSr' DAMC01*E2 
2-PI * RHO * E3) 

ClKl = E3 



RETURN 
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C THE RK21 SUBROUTINE DETERMINES THE K3 COEFFICIENTS FOR THE 

C COMBINED MASSES 

SUBROUTINE RK21 

COMMON Rl(512) ,R2(512) ,C1 (512) ,C2(512) ,TAU(512) , PI , AMU,GAMMA1 ,PSI , 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R1K1,R2K2,R2K3, 
3C2K0,C2Kl,C2K2,C2K3,L,DAMCOl,DELTAU,F01 
B = DELTAU / 2. 

T = TAU(L) + B 
FI = 2. * PI * GAMMAl * T 
El = R1(L) + B * RlKl 
E2 = R2(L) +B * R2K1 
E3 = C1(L) + B * cm 
R1K2 = E2 

R2K2 =(4.*PI/(1.+AMU))*(PI*SIN(FI)-DAMC01*E2-PSI*DAMC01*E2-PI*E1- 
2PI* RHO * E3) 

C1K2 = E3 
RETURN 



END 
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THL KOI SUEROUTINi: DETERMINES THE K4 COEFFICIENTS FOR THE 
COMBINED MASSES 
SUBROUTINE RK31 

COMMONRl(512) ,R2 (512) ,C1(512) ,C2(512) ,TAU(512) ,PI,AMU,GAMMA1,PSI, 
2RHO,R1KO,R1K1,R1K2,R1K3,C1KO,C1K1,C1K2,C1K3,R2KO,R2K1,R2K2,R2K3, 
3C2KO,C2K1,C2K2,C2K3,L,DAMC01,DELTAU,F01 
B = DELTAU 
T = TAU(L) + B 
FI = 2. * PI * GAMMAl * T 

El = R1(L) + B * R1K2 

E2 = R2(L) + B * R2K2 

E3 = C1(L) + B * C1K2 

R1K3 = E2 

R2K3 =C4.*PI/(l.+AMU))*(PI*SIN(FI)-DAMC01*E2-PSI*DAMe01 *E2-PI* 
2RH0 * E3) 

C1K3 = E3 



RETURN 
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C THE AREA FUNCTION COMPUTES: Y(N+1) - Y(N) FOR R1,R2,C1,C2 

FUNCTION AREA(B1,B2,B3,B4,B5,B6) 

AREA = B1 + B2/6. * (B3 + 2.*B4 + 2.*B5 + B6) 

RETURN 

END 
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APPENDIX B 



In this appendix the individual platform system responses 
are presented. Each plot is serialized below the abscissa 
in a statement, • Run Code XYZ." The code is interpretted 
as follows; X = 1; the plot is the result of a R-K integra- 
tion. 

Y = 1 through 7; denotes the value of Yi used in 
the R-K integration. 

Z = 1 through 3; denotes the value of FOl used 
during the R-K integration. 

A complete description of the codes is listed in Table B.l. 

TABLE B.l, Computer Run Code Summary 



CODE 


Yi 


FOl 


CODE 




FOl 


111 


0.5 


0.01 


152 


2.0 


0.001 


121 


0 . 8 


0.01 


162 


3.0 


0. 001 


131 


1 . 0 


0.01 


172 


5 . 0 


0.001 


141 


1 . 5 


0.01 


113 


0 . 5 


0 . 0001 


151 


2.0 


0.01 


123 


0 . 8 


0.0001 


161 


3.0 


0.01 


133 


1.0 


0.0001 


171 


5 . 0 


0.01 


143 


1.5 


0.0001 


112 


0.5 


0.001 


153 


2.0 


0.0001 


122 


0.8 


0.001 


163 


3.0 


0. 0001 


132 


1 . 0 


0 . 001 


173 


5 . o' 


0.0001 


142 


2.0 


0.001 
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Each plot is followed by an analysis of the platform response 
with respect to steady-state magnitude, the effect of the 
excitation frequency, the position of the dynamic equilibrium, 
positive or negative, relative to the static equilibrium zero 
position and the primary response frequencies as determined by 
a Fourier transform of the platform response. The plots which 
were computed for FOl equal to 0.0001 are included without 
an explanation of the plots. These plots coincided with the 
platform response plots for FOl equal to 0.001 for all values 
of Yi* Therefore no additional information could be obtained 



from them. 
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FIGURE B.l, Platform Response for FOl = 0.01; Yi* 0-5 
The steady-state platform displacement is 3.5 units. 

The period between maximum displacements is 2x which corres- 
ponds to the excitation frequency of 1 radian/second. The 
dynamic equilibrium position is slightly positive. The Fourier 
transform of the platform response yields primary frequencies 
of 0.42 and 0.84 radians/second. These frequencies correspond 
to the difference between the excitation frequency and the 
conductor natural frequency, and the first harmonic of this 
frequency . 
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FIGURE B.2, Platform Response for F01= 0.01; Yi = 

The steady-state magnitude of the platform response is 
8.87 units. The period of the response is 1.25x corresponding 
to the excitation frequency. The dynamic equilibrium position 
is 3.84 units negative. The primary frequencies of the re- 
sponse are the excitation frequency, the sura of the excitation 
frequency and the platform natural frequency, and an unident- 
ified combination of frequencies equal to 5.7 rad ians/ second . 
More comprehensive data is required to determine the exact 
components of this frequency. 
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TXIC PARMCTER RUN CODE 131 



FIGURE B.3, Platform Response for F01= 0.01; Yi= 1.0 

The average steady-state response magnitude is 5.19 units. 
The frequency of the displacement is equal to the excitation 
frequency. As the overall displacement moves in the negative 
direction, the individual cycle amplitudes increase over a 
four cycle period such that if the equilibrium position were 
normalized, a beat frequency would be prevalent. If this 
conclusion is true, the period over the four cycles should be 
equal to one-half the beat frequency period, or . 

Substituting the natural frequencies of the platform and 
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the conductor for Wj and Uj respectively, the period is 4.28t. 
By inspection of the displacement plot, the period of the 
four cycles is approximately 4 t. Thus, a beat frequency 
caused by the interaction of the platform and the conductor 
is a plausible explanation of the platform displacment plot. 
The dynamic equilibrium position is 2.36 units negative. The 
Fourier transform of the response indicates a primary 
frequency equal to the excitation frequency and the platform 
natural frequency. 
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FIGURE B.4, Platform Response for F01= 0.01; y y- 1.5 
The steady-state response magnitude is 0.67 units. The 
period over one response cycle is 0.66t corresponding to the 
excitation frequency. However, the period between two adjacent 
peaks of equal amplitude is 2t. The dynamic equilibrium 
position is slightly positive. From the pourier transform of 
the response, the principle frequencies are the excitation 
frequency, the platform natural frequency and a frequency equal 
to the difference between these two frequencies. The period 
of the third frequency is 2t corresponding to the predominant 
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period between displacement peaks. 
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FIGURE B,5, Platform Response for F01= 0,01; 2.0 

The steady-state response magnitude is 0,50 units. The 
platform oscillates alternately with a period which equates to 
the sum of the platform and conductor natural frequencies and 
a period which is twice the summation frequency. The variation 
in response period implies that the platform and conductor are 
alternately in and out-of-phase. The dynamic is slightly 
positive. The Fourier transform indicates primary frequencies 
equal to the excitation frequency and the platform natural 
frequency . 
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FIGURE B.6, Platform Response for F01= 0.01; y^= 3.0 
The average steady-state magnitude is 0.2 units. In 
steady-state, the platform response has a beat frequency 
superimposed. The beat period is 20x which is a frequency 
of 0.1 radians /s econd . The dynamic equilibrium position is 
slightly positive. The primary response frequencies are the 
excitation frequency and the platform natural frequency. 
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FIGURE B.7, Platform Response for F01= 0.01; Yi= 5.0 
The steady-state magnitude is 0.012 units. Uniquely, 
thf system settled into an unstable equilibrium position 
before reaching its final, steady-state equilibrium position. 
The system cycles at a frequency of 2.65 radians/second 
which equates approximately to the sum of the platform and 
the conductor natural frequencies. It appears that at this 
high excitation frequency the relative velocities of the 
conductor are at a level such that the platform and conductor 
momenta are nearly equal. When the conductor and platform 
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are out -o f -phase , the momentum of the conductor is sufficient 
enough to cancel the momentum of the platform causing a 
reduction in platform motion. When the two components are 
in phase, the platform displacement nearly doubles. The 
dynamic equilibrium position coincides with the static 
equilibrium position. The pourier transform produces frequen- 
cies from the platform response equal to the conductor natural 
frequency, a frequency which is twice the conductor natural 
frequency and a frequency equal to the sum of the platform 
and the conductor natural frequencies. 
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FIGURE B.8, Platform Response for F01= 0.001; Yi= 0*5 
The steady- state magnitude is 1.35 units. The period 
of the response is 2x which equates to the excitation 
frequency. The dynamic equilibrium position is slightly 
negative. From the Fourier transform, the primary response 
frequencies are the excitation frequency and the platform 
natural frequency. 
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FIGURE B.9, Platform Response for F01= 0.001; Yi= 0-8 
The average steady-state magnitude is 2,6 units. The 
individual displacement cycles have a period of 1.25 t which 
is the excitation frequency. There is a slight beat frequency 
superimposed which has a period of IOt and a frequency of 
0.2 radians/second. The dynamic equilibrium position is 
slightly positive. The Fourier transform yields primary 
frequencies equal to the excitation frequency and the platform 
natural frequency from the platform response. 




FIGURE B.IO, Platform Response for F01= 0.001; Yi= 1*0 
The average steady-state magnitude is 5.3 units. The 
platform displacement period is It corresponding to the 
excitation frequency and the platform natural frequency. 

A beat frequency is more porminent . The beat period i.s 4x 
wh:;.ch is a frequency of 0.52 radians/ second . The dynamic 
eqtiilibrium position is 0.55 units negative. From the 
Fourier transform, the primary frequency is the excitation 
frequency. Nothing in the transform indicates the component 
frequencies which make up the beat frequency. 
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FIGURE B.ll, Platform Response for F01= 0.001; Yi= I--’ 

The average steady-state magnitude is 1.1 units. ’he 
period of a displacement cycle is 0.67 t which equates to the 
excitation frequency. The beat phenomena has a period of 
24 t and a frequency of 0.08 radians/second . The dynamic 
equilibrium position is 0.07 units positive. The Fourier 
transform of the platform response indicates primary frequen- 
cies equal to the excitation frequency, the platform natural 
frequency and a frequency of 1.87 radians/second which requires 
further analysis for component identification. 
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FIGURE B.12, Platform Response for F01= 0.001; y^= 2.0 
The steady- state displacement magnitude is 0.50 units. 

The period of the response cycle cannot be uniquely determined 
as the period varies over adjacent cycles as the displacement 
amplitude alternately increases and decreases. However, over 
four cycles the average period is 0.5 t which corresponds to the 
excitation frequency. The period of the amplitude alteration 
is 4 t which is a frequency of 0.53 radians/second. The 
dynamic equilibrium position is negative slightly. The 
Fourier transform of the response indicates primary frequen- 
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cies equal to the excitation frequency and the platform 
natural frequency. 
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FIGURE B.13, Platform Response for F01= 0.001; Yi“ 3.0 
The average steady-state magnitude is 0.2 units. The 
displacement period is 0.33 t which is the excitation frequency. 
The platform response exhibits a beat phenomena with a period 
of 24 t and a frequency of 0.08 radians/second . The dynamic 
equilibrium position is slightly positive. From the Fourier 
transform of the response, the primary frequencies are the 
excitation frequency and the platform natural frequency. It 
appears that the excitation frequency is sufficiently removed 
from the platform natural frequency that excitation force 
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is an insignificant factor. This response plot duplicates 
the response plot for FOl equal to 0,01. 
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FIGURE B,14, Platform Response for F01= 0.001; Yi= 5,0 
This platform response plot coincides exactly with 
the plot for FOl equal to 0.01. The steady-state magnitude 
is 0.012 units. The dynamic equilibrium position coincides 
with the static equilibrium position. The primary frequen- 
cies from the Fourier transform of the response are the 
sum of the platform and conductor natural frequencies, the 
conductor natural frequency and a frequency which is twice 
the conductor natural frequency. 
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FIGURE B.15 



Platform Response for F01= 0.0001; Yi= 0.5 
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B.16, Platform Response for F01= 0.0001; Yjs 0*8 
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FIGURE B.17, Platform Response for F01= 0.0001; Yj= 1-0 
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FIGURE B.18, 



Platform Response for F01= 0.0001; y 1.5 
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FIGURE B. 19, Platform Response for F01= 0.0001; y^= 2.0 
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B.20, Platform Response for F01= 0.0001; y = 
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